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Abstract 

We investigate the signals of a supersymmetric scenario where the lighter tau- 
sneutrino is the lightest supersymmetric particle, while the lighter stau-state is the next 
lightest. We confirm that such a scenario can be motivated within the framework of 
minimal supergravity, with just the addition of a right-chiral neutrino superfield. Such 
a spectrum leads to rather unusual signals of supersymmetry, showing stable tracks 
of the stau in the muon chambers. We have studied two types of signals, namely (a) 
two or more hard jets and two stau tracks, and (b) hard jets with two muons and two 
stau tracks. We demonstrate that the stau tracks can be distinguished from the muonic 
ones through proper kinematic cuts which also enable one to remove all standard model 
backgrounds. 
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1 Introduction 



Supersymmetry (SUSY) is still the most sought-after new physics option around the TeV 
scale, and concerted efforts are on to discover it in any of its possible ramifications at the 
Large Hadron Collider (LHC). Apart from stabilizing the electroweak symmetry breaking 
sector and providing a rather tantalizing hint of Grand Unification, SUSY (in the R-parity 
conserving version) also provides a cold dark matter candidate in the form of the stable 
lightest supersymmetric particle (LSP). The most common practice is to assume that the 
lightest neutralino is the LSP [1], and the 'standard' signals such as (jets + fir), (dileptons 
+ fix) etc. are widely studied under this assumption [2]. After decades of rather exhaustive 
investigation, two questions that are now being asked are, (a) how tractable are SUSY signals 
with other types of LSP? And (b) is it possible to have quasi-stable but cosmologically 
allowed charged particles in a SUSY scenario, with distinct collider signatures? As we shall 
see below, the answers to both questions can be of diverse nature but nonetheless related, 
and are worth exploring in detail before the LHC takes off. 

It must be realized at the same time that hints of new physics stare us in the face at a 
scale much lower than a TeV, as one finds incontrovertible evidence of neutrino masses and 
mixing of rather unexpected nature [3]. The shortest step beyond the standard electroweak 
model, which provides an explanation of them, consists in hypothesizing a right-handed 
neutrino for each fermion family [4]. If one accepts such a hypothesis, and at the same 
time depends on SUSY for the stabilization of the electroweak scale, then three additional 
members are added to the assortment of R-odd (super)particles, namely, the right-chiral 
sneutrinos. Interestingly, such sneutrinos are perfectly inactive with respect to all the gauge 
interactions, and have only couplings proportional to the neutrino masses. Of course, the 
physical states in each family involve a tiny mixing of both chiral types. The non-interactive 
nature of the dominantly right-chiral state relative to the left-chiral one makes it possible 
for the former (as opposed to the latter [5]) to be a dark matter candidate [6]. Therefore, it 
makes perfect sense to envision a SUSY spectrum with such a sneutrino as the LSP. 

In this paper, we study some signals of such a dominantly right-sneutrino LSP (ui) at 
the LHC. When this happens, a number of candidates exist for the next-to-lightest SUSY 
particle (NLSP). These include the dominantly left-sneutrino state (z> 2 ), a slepton (especially 
the stau, with its opportunity to have a rather light mass eigenstate with large tan (3, tan (5 
being the ratio of the vacuum expectation values of the two Higgs doublets), a neutralino, 
a chargino, or even a stop in special cases. Among these, we consider the particular case 
of a stau NLSP because the collider signals are rather typical in the pronounced difference 
with the neutralino LSP scenario. The consequences of a spectrum with, for example, a stop 
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NLSP has been studied in the literature [7]. 

Interestingly, for a v x LSP, any NLSP will decay into it only through an interaction 
strength proportional to the neutrino mass. This is because Yukawa interactions will invari- 
ably involve this mass, while gauge interaction will depend on the admixture of the SU(2) L 
doublet component in if\ proportional to the neutrino Yukawa coupling. Thus the decay of 
the NLSP to the LSP is too suppressed over most of the viable parameter space to take place 
within the detector. While the LSP still contributes (at least partially) to cold dark matter, 
the NLSP appears to be stable, as far as collider detectors are concerned. Consequently, 
SUSY signals in the case considered by us will be marked by charged tracks characteristic 
of staus. 

'Stable stau' phenomenology of the kind mentioned above can in general be of several 
origins, namely, 

• Gravitino LSP in a supergravity (SUGRA) theory [8]. Here the gravitino is of such 
mass as to be a cold dark matter candidate, and its coupling to the stau (being inversely 
proportional to the gravitino mass) is suppressed enough to prevent the latter from 
decaying within the detector. 

• Stau-NLSP in gauge mediated SUSY breaking (GMSB) models [9]. Here the stau may 
be stable if the effective F-term leading to the dynamical generation of gravitino mass 
is on the higher side (above about (10 7 GeU) 2 ). 

• A neutralino LSP nearly degenerate with a stau NLSP [10]. This is the so-called co- 
annihilation region, with the small mass difference suppressing the decay of the NLSP. 

• An LSP dominated by the right-chiral sneutrino [11], with a stau NLSP. As we shall 
see, such a possibility is quite viable in a SUGRA scenario admitting a right-neutrino 
superfield in the observable sector. 

The present study concerns the last one of the above scenarios. Although long-lived stau 
NLSP's arising from one or the other of the possibilities mentioned above are frequently 
discussed nowadays, our purpose is to specially emphasize the following points: 

(a) The stau can become an NLSP, and the corresponding right sneutrino, the LSP, in 
rather natural regions of the SUGRA parameter space, provided that right sneutrino mass 
is allowed to evolve from the common scalar mass at high scale. 

(b) The stau NLSP is likely to leave charged tracks on reaching the muon detector [12]. 
Ways of distinguishing it from a muon have been suggested on the basis of the thickness 
of the track, and also through the time delay between the inner tracking chamber and the 
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muon detector. Moreover, the use of stoppers has been sometimes advocated for intercepting 
long-lived staus, whereby their late decays can be studied in underground chambers at a later 
stage [13]. While the efficacy of such methods is admitted, we wish to point out that the very 
kinematic properties of the tracks in the muon chamber set the long-lived staus apart from 
the muons in a conspicuous fashion, and the characteristic signal events of such a scenario 
can be separated from the standard model (SM) backgrounds in a rather straightforward 
manner. 

In section 2 we discuss how the z>i-LSP, stau-NLSP scenario can arise naturally within the 
SUGRA framework including a right-chiral neutrino superfield for each generation. Section 
3 contains discussions on the two principal signals of such a scenario, which can originate 
in the production and cascade decays of strongly interacting superparticles. These are the 
(hard jets + f-pair) and (hard jets + f-pair + dimuon) signals, respectively. 1 . We summarize 
and conclude in section 4. 



2 Right sneutrino LSP in supergravity 

The simplest extension to the SM spectrum, which would explain non-vanishing neutrino 
masses [14], is to add one right-handed neutrino for each generation. However, having such 
small Dirac masses for the neutrinos would imply that the neutrino Yukawa couplings are 
quite small. Considering the different neutrino mass hierarchies as benchmark scenarios, a 
normal (or inverted) hierarchy would imply that such couplings are at most of the order of 
~ 1CT 13 , whereas they can be larger by an order or so for the case of degenerate neutrinos. 

Here we consider a scenario where lepton number is conserved. Thus the superpotential 
of the minimal SUSY standard model (MSSM) is extended by just one term which, for a 
particular family, is of the form 

W* = y v H u U c R (1) 

where y v is the Yukawa coupling, L is the left-handed lepton superfield and H u is the Higgs 
superfield responsible for giving mass to the T3 = +1/2 fermions. The above term in the 
superpotential obviously implies the inclusion of right-handed sneutrinos in the particle spec- 
trum. As discussed in the introduction, such sneutrinos will have all their interactions pro- 
portional to the corresponding neutrino masses. So the dominantly right-handed eigenstate 

1 Although di-electron signals are as important as dimuons, we demonstrate our predictions with reference 
to the latter only, since the stable staus may appear like muons. Hence the challenge of separating such 
'fake' muons with real ones is best reflected in a analysis involving dimuon samples. 
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of the tau-sneutrino might become a possible candidate for the LSP. Working in the frame- 
work of minimal supergravity (mSUGRA) model of SUSY, we find that one can have the 
lightest sneutrino (composed mostly of the right-chiral component) as the LSP, consistent 
with all experimental bounds [16] and also within the acceptable limits of dark matter density 
in the universe [11, 17]. 

The neutrinos obtain mass as: 



The actual mass eigenvalues will of course depend on the Yukawa couplings which may 
have a hierarchy among them, and will in general be matrix-valued. However, the different 
neutrino mass hierarchies, and the consequent hierarchies in the Yukawa couplings, do not 
affect the generic feature of the collider signal under consideration here, as the decay rate 
of the stau-NLSP is extremely suppressed in either of them. It is basically the smallness of 
the Yukawa coupling which plays a crucial role in our understanding of the spectrum and 
its consequent features. 

Upon inclusion of right-chiral neutrino superfield into the SUGRA fold, the superparticle 
spectrum mimics the mSUGRA spectrum in all details except for the identity of the LSP. 
SUSY breaking effects are introduced explicitly as universal soft terms for scalars (m ) and 
gauginos (^1/2) together with the so-called A and B-parameters (of which the latter is 
determined by electroweak symmetry breaking conditions) in the Lagrangian. Masses of 
squarks, sleptons and gauginos, all the mass parameters in the Higgs sector as well as the 
Higgsino mass parameter \x (upto a sign) are determined, once the SUSY breaking parameters 
at a high scale (O ~ 10 11 GeV) are specified. The mass terms for sneutrinos (neglecting 
inter-family mixing) arising in this manner are given by 



where A v is the term driving left-right mixing in the scalar mass matrix, and is obtained 
by running of the trilinear soft SUSY breaking term A. The Yukawa couplings can cause 
large splitting in the third-family squark and sleptons masses while the first two families 
are effectively degenerate. On the other hand, one expects minimal left-right mixing of 
sneutrinos as the Yukawa couplings are all extremely small. 
The mass-squared matrix for the sneutrino thus looks like 



where Mi is the soft scalar mass for the left-handed sleptons whereas the M„ R is that for the 
right-handed sneutrino. In general, 7^ M i/R because of their different evolution patterns 




(2) 



£ sof t ~ M? R \i> R \ 2 + (y u A v H u .LI> c R + h.c.) 



(3) 
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as well as the D-term contribution for the former. While the evolution of all parameters of 
minimal SUSY remain practically unaffected in this scenario, the right-chiral sneutrino mass 
parameter evolves [6] at the one-loop level as: 

dt ~ IGtt 2 ^ Av ■ [b) 

Clearly, the extremely small Yukawa couplings cause M„ R to remain nearly frozen at the 
value mo, whereas the other sfermion masses are jacked up at the electroweak scale. Thus, 
for a wide range of values of the gaugino mass, one naturally has sneutrino LSP's, which, 
for every family, is dominated by the right-chiral state: 

z>i = — v L sin 9 + u R cos 9 (6) 

The mixing angle 9 is given as 

tan2g= 2 ^ Sin f |CQt/ ! /i -^ 1 (7) 

which is clearly suppressed due to y v . It is to be noted that all three (dominantly) right 
sneutrinos have similar fate here, and one has near-degeneracy of three such LSP's. However, 
of the three charged slepton families, the amount of left-right mixing is always the largest in 
the third (being, of course, more pronounced for large tan/3), and the lighter stau (fx) often 
turns out to be the NLSP in such a scenario. 2 

Thus the mSUGRA parameter set (mo, mi/ 2 , A, sign(/i) and tan/3) in an R-parity con- 
serving scenario can eminently lead to a spectrum where all the three generation of right- 
sneutrinos will be either stable or metastable but very long-lived, and can lead to different 
decay chains of supersymmetric particles, as compared to those with a neutralino LSP. How- 
ever, as we shall see below, the controlling agent is the the lighter sneutrino eigenstate of 
the third family, so long as the state f x is the lightest among the charged sleptons. 

We have chosen the above scenario to demonstrate our main point for another reason. 
In general, apart from the stau, either a neutralino or a chargino (being the remaining R- 
odd non-strongly interacting particles) could be the NLSP. However, a neutralino NLSP will 
decay into a neutrino and the sneutrino LSP, both of which are invisible. Thus the collider 
signals will not be qualitatively different from those of a neutralino LSP model, irrespective 
of whether the NLSP decays inside the detector or outside. A chargino NLSP, on the other 



2 We have neglected inter-family mixing in the sneutrino sector in this study. While it is true that 
near-degenerate physical states makes their mixing a likely consequence, the extent of such mixing is model- 
dependent, and does not generally affect the fact that all cascades culminate in the lighter stau, so long as 
the latter in the NLSP, which is the scenario under study here. 
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Figure 1: Lifetime of stau NLSP against the universal gaugino mass parameter m 1( / 2 . Solid (red) 
line is for tan/3 = 10 and the dashed (black) line is for tan(3 = 30. Other SUGRA input parameters 
are: m = 100 GeV, A = 100 GeV, sgn(fi) = 1. 

hand, can make a difference in this respect through charged tracks. However, it is difficult to 
construct a model where the lighter chargino (xt) is lighter than the lightest neutralino (xi) 
unless gaugino non-universality of a rather gross nature is postulated. We have chosen to 
confine ourselves to a relatively minimal high-scale picture, where such an option is somewhat 
far-fetched. The remaining situation is one with a f NLSP, which is both (as we have shown 
above) easy to accommodate within an mSUGRA scenario and gives qualitatively different 
SUSY signals. 

Any (super) particle will have to decay into the lighter sneutrino of a particular family via 
either gauge interactions (such as II — > Wv\) or Yukawa coupling (such as II — > H~V\ or 
z>2 — ► hPvi). In the former case, the decay depends entirely on the vl content of i/i, which 
again depends on the neutrino Yukawa coupling. The same parameter explicitly controls the 
decay in the latter case, too. Therefore, while the lighter sneutrinos of the first two families 
can in principle be produced from decays of the corresponding charged sleptons, such decays 
will be always suppressed compared to even three-body decays such as ei(fli) — > e(/i)ffi. 
For the f-NLSP, however, the only available decay channel is f\ — > Wui, which is driven 
by the extremely small neutrino Yukawa coupling. 

This causes the NLSP to be a long-lived particle and opens up a whole set of new 
possibilities for collider signatures for such long-lived particles, while retaining contributions 
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to cold dark matter from the sneutrino LSP 3 . The NLSP appears stable in collider detectors 
and gives highly ionizing charged tracks. 

With all the above considerations in mind, we concentrate on the lighter stau (fi) to be the 
NLSP with lifetime large enough to penetrate collider detectors like the muons themselves. 
Using the spectrum generator of the package IS A JET 7.69 [18], we find that a large mSUGRA 
parameter space can realize this scenario of a right-sneutrino LSP and stau NLSP, provided 
that m < my 2 an d one has tanfi of the order of 10 and above, the latter condition being 
responsible for a larger left-right off-diagonal term in the stau mass matrix (and thus one 
smaller eigenvalue). 

It is worth mentioning that the region of the mSUGRA parameter space where we work is 
consistent with all experimental bounds, including both collider and low-energy constraints 
(such as the LEP and Tevatron constraints on the masses of Higgs, gluinos, charginos and 
so on as well as those from b — > 57, correction to the p-parameter, (g^ — 2) etc.). Our 
choice of parameters in the mo — mi/2 plane would correspond to a stau-LSP without the 
right-sneutrino in the (super)particle spectrum, which is ruled out in the context of x?-LSP 
scenario [19]. However with the right-sneutrino as the LSP, we find this choice to be a 
preferable and well-motivated option. The z>i-LSP arising out of such a choice becomes a 
viable constituent of cold dark matter, though not necessarily the only one. 

We mostly focus on the region of the parameter space where m^ > m^ 1 + mw is satisfied, 
so the dominant decay mode is the two-body decay of the NLSP, fi — > v{W . In Figure 1 
we show the lifetime of such long-lived stau's plotted against the universal gaugino mass 
parameter m 1 / 2 for a particular choice of (m , A, sign(fi) and tan/3). One must however 
note that having too large a value for my 2 will make the gluinos and squarks too heavy 
to be produced at colliders. We have chosen the case of degenerate neutrinos, since this 
permits the largest possible values of the Yukawa couplings. The fact that the lifetimes are 
so large even with this choice is a convincing demonstration of the long-lived nature of the 
stau-NLSP, as far as collider detectors are concerned. In Table 1 we present two benchmark 
points for our study of such long-lived staus at the LHC. In the next section we discuss 
the signatures of the stau NLSP at the LHC, concentrating on the two types of final states 
mentioned in the introduction. Using the formulae given by Moroi et al. in reference 4, 
the contribution to the cold dark matter relic density is found to be about one order of 
magnitude below the acceptable value of Vth 2 . While this leaves room for additional sources 
of dark matter, the scenario presented here is consistent from the viewpoint of over-closure 
of the universe. 

3 In fact, the possibility of accounting at least partially for the cold dark matter of the universe makes the 
right-sneutrino LSP scenario slightly more viable than one with a supcrlight gravitino playing a similar role. 
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3 Signatures of stau-NLSP at LHC 



In this section we discuss the signatures of the long-lived stau-NLSP at the LHC and con- 
centrate on two different final states, viz. 

• 2fi + 2(or more) jets (p T > 1 00 GeV) 

• 2fi + dimuon + 2(or more) jets (px > 100 GeV) 

Keeping the above signals in mind, we focus on the two benchmark points listed in Table 1 
and study their signatures at the LHC, for an integrated luminosity of 30 fb~ l . 
The low energy mass spectrum schematically takes the following form: 

< m fl < m^o < < ... < m~ g 

which would suggest that all superparticle productions at the LHC would finally end at the 
sneutrino LSP. However, as discussed in the previous section and also evident from Figure 1, 
the lifetime of the fi-NLSP is quite large, ranging between a few seconds to a few minutes, 
making it stable in the context of collider studies. So it will almost always decay outside the 
detector, leaving characteristic signals like charged tracks, with large transverse momenta. 
In fact this would be quite a contrast to the traditionally thought of SUSY signals with 
large missing transverse energy. This fact acts as the main thrust of our work, as the stable 
stau will behave just like a muon. However, in the absence of spin identification, these 
staus will behave more like very heavy muon-like particles with v/c) < 1 and such 
heavy charged particles will have high specific ionisation due to their slow motion within the 
detector. Studies exist in the literature which have looked at such heavy and stable charged 
particles and considered their ionisation properties and time of flight as distinctive qualities 
which separate them from the muons [12]. In this work we take a qualitatively different 
approach towards studying the long-lived staus which would be based more on an analysis 
pertaining to studying the kinematics of processes producing such particles at the LHC. We 
will see through our analysis, that by simply looking at certain kinematic distributions, one 
can distinguish such long-lived staus from the muons. In the kinematic plots presented later 
in this section, we have neglected the energy loss of the staus reaching the muon chambers. 
Including this loss is matter of detailed detector analysis, and are unlikely to change the 
general features of the tracks, which are qualitatively rather distinctive. 

It should be noted that small variations of the values of parameters used here may lead 
to one selectron and smuon state each to be also lighter than the lightest neutralino. In 
that case, although the only two-body decay available to each of them lead to the lighter 
sneutrino states of the first two families, such decays are tremendously suppressed due the 



9 



Parameter 


Rpnchmark noint 1 

J 9 V.' X X V ' l X l X X CXI XV Y-* X X X Li _l_ 


Rpnchmark noint 2 








m n = 100 GeV m, /9 = 600 GeV 


mn = 110 GeV m, /9 = 700 GeV 


mSUGRA input 


A = 100 GeV, s#n(/i) = + 


A = 100 GeV, s#n(/i) = + 








tan /? = 30 


tan /3 = 10 




// 




694 


810 


m e~ L i m fi L 


420 


486 


m e~R ? m fi.R 


9^1 


9SQ 


m i'eL J m i' IJ ,L 


419 


47Q 






478 




100 


110 


m fl 


187 


281 


TYlf 2 


422 


486 


Xi 


9/n 




m v o 


4fiQ 






i uu 


81 ^ 

Oil) 




713 


829 




470 


552 


m xt 


713 


829 


m-g 


1366 


1574 




1 9^7 


1494 




1 1 


1 ^7^ 




1 9^Q 








1 3fi7 


m h 


Q84 


1 1 37 


m h 


1176 


1365 


m h 


1123 


1330 


m b 2 


1161 


1358 


m h o 


1 18 


1 18 


m H o 


712 


941 


rriAo 


707 


935 


m H ± 


717 


944 



Table 1: Proposed benchmark points for study of stau-NLSP scenario in the SUGRA fold with right- 
sneutrino LSP. All superparticle masses are given in GeV. Due to very small mixing, i>\ — vr. The 
top mass, mt = 111 A GeV [20] has been used for running the parameters. 
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same reason as in the case of the stau-NLSP. On the other hand, each of them can have 
three-body decay, mediated, for example, by a virtual neutralino, into a stau, a tau and 
an electron(muon). This channel will prevent the selectrons(smuons) from being long-lived, 
while they will give rise to additional events with stau charged tracks. 

We now take up the above mentioned signals one by one and present our results in the 
following subsections. 

3.1 2fi + hard jets 

For our simulation we have generated 10 million events for SUSY particle production, and 
allowed their decay through all possible modes in to lighter stau. The signals mentioned here 
arise mostly from the direct decay of gluinos and squarks, produced via strong interaction, 
into the lightest neutralino, with the latter decaying into a tau and a lighter stau, and 
the tau decaying hadronically in turn. Cascades through other neutralinos and charginos 
supplement the rates to a moderate extent. We have used PYTHIA 6.409 [21] for our 
event generation and interfaced it with ISASUGRA, the SUSY-spectrum generator routine 
contained in ISA JET 7.69. Following the procedure outlined in the previous section, we have 
introduced the right-sneutrino in the ISAJET program to generate the mSUGRA spectrum 
with a right-sneutrino LSP with its corresponding renormalization group equation RGE (eqn. 
5) also included. To make the Xi decay in PYTHIA we have modified the decay table to 
include its decay to f\f and demand that the Xi is n °t the LSP. 

The parton densities have been evaluated at Q — 2m^ using CTEQ5L [22], and the 
renormalization scale and factorization scale are 

= Q = Hr 

(On setting the scale at the squark/gluino masses the rates are reduced by about 20%, and 
none of the general conclusions is altered). The effects of both initial state radiation (ISR) 
and final state radiation (FSR) have also been taken into account. We have included the 
effects of hadronization and multiple interaction with the help of PYTHIA hadronization 
schemes. 

To define jets we use the simple-minded jet cone algorithm implemented in PYTHIA 
through the subroutine PYCELL. The basic parameters we have used in this cluster-finding 
routine are as follows: 

• The jet conical width is ARjj > 0.7, where 

AR n = yjA^ + 

where Arjjj and Acpjj are intervals in rapidity and azimuthal angle. 
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• The Gaussian smearing of the E± with a width proportional to 1.2E± is used to take 
into account the energy resolution of the detector. 

• The summed E± of the jet (consisting of all cells within the cone of radius R in the 
(rj, 4>) plane) should be greater than 50 GeV to be accepted as a jet. 

• The rj coverage range for jets is taken to be \rj\ < 3.0. 

To select our final states, we demand the following requirements (called basic cuts) on our 
sample events: 

• Each f\ should have pt > 30 GeV. 

• Both the fi's should satisfy \r]\ < 2.5, to ensure that they lie within the coverage of 
the muon detector. 

• ARf lfl > 0.2, to ensure that the fi's are well resolved in space. 

• At least two jets with p T > 100 GeV (hard jets). 

• In addition we have rejected events having photons with |?y 7 | < 2.5 and > 25 GeV. 
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Figure 2: Kinematic distributions for the signal 2ti + (> 2) hard jets. In (a) the transverse momen- 
tum distributions for the harder f\ is shown and (b) shows the invariant mass distribution for the f\ 
•pair. The dash-dot-dash (red) histograms are for benchmark point 1 and the solid (blue) histogram 
for benchmark point 2. The dashed (black) histograms show the corresponding SM background. 
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Figure 3: Kinematic distributions for the signal 1f\ + (> 2) hard jets after imposing the stronger 
cut, pt > 350 GeV on both charged tracks. In (a) the transverse momentum distributions for the 
harder f\ is shown and (b) shows the invariant mass distribution for the t\ pair. We follow the 
same notation as in Figure 2. 

As the charged tracks pass through the muon chamber, it is probable that the muonic events 
will fake our signal. Therefore, events with two or more hard jets and two central muons 
will prima facie constitute our standard model background. The leading contribution to 
such final states satisfying our basic cuts comes from top-pair production and its subsequent 
decay into dimuons, with similar topology as that of the signal. The sub-leading contributions 
consist in weak boson pair production. For performing the background analysis we again 
used the same criteria as stated above for the signal, with all kinematic features of the 
long-lived charged tracks attributed to the muons. 

We would like to emphasize here that the demand for two (at least) hard jets with 
Pt > 100 GeV may itself be expected to act as a 'killer' to background, as we are well 
aware of the fact that events with large number of high p? jets is itself a 'hint' to SUSY 
provided that we do not have other candidates for physics beyond the SM. However, SUSY 
signals have this advantage ordinarily through the occurrence of large missing- /p?. In this 
case, the hard stau tracks assume the role of the latter. The other important point to 
note is that the background is almost completely reducible with the imposition of stronger 
event selection criteria, as we shall see below. In Figure 2, we present distributions of a 
few observables where we could distinguish the signal from backgrounds. These are the 
transverse momentum (pt) of the harder charged track (and the harder muon in the case 
of backgrounds) and the invariant mass (M 12 ) of the two charged tracks (or dimuons). In 
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addition, the radii of curvatures of the stau tracks will also lie in a clearly distinguishable 
range, as seen from the p^-distributions. 



Cuts 


Background 


Benchmark point 1 


Benchmark point 2 


Basic 


39617 


8337 


1278 


Basic + p T > 350 GeV 


5 


2587 


737 



Table 2: The expected number of events for the signal and background with the different cuts 
imposed on the selection of events. We have assumed an integrated luminosity of 30 fb~ l at the 
LHC to generate the events. 



Figure 2 indicates that the basic cuts retain a large number of background events. A 
substantial fraction of these events come from the Z-peak, as evident from the invariant 
mass distributions. However, on examining the p^-distributions, we find it most convenient 
to eliminate them by imposing a stronger p? cut on both the charged tracks. This cut has 
been set at 350 GeV. The effect is rather dramatic, and makes the signal stand out clearly 
for both the benchmark points, as can be seen from Table 2. 

Figure 3 shows the corresponding distributions after the backgrounds are nearly gone. 
These constitute the real 'signal distributions', and can be used to extract information, for 
example, about the mass of the stau-NLSP and other SUSY parameters. 

Finally, one has to remember that, in addition to the removal of SM backgrounds, one 
needs to distinguish the signal mentioned above from that of a SUSY scenario where the 
lightest neutralino is the LSP. There, too, a pair of high-p^ muons can arise from cascades, 
together with hard jets. However, as has already been mentioned, the neutralino LSP causes 
the missing-p T distribution to be much harder, at least when the LSP mass is related through 
universality conditions to the higher charginos and neutralinos which, in turn, are constrained 
by results from the Large Electron Positron (LEP) collider. Consequently, signals in the 
scenario under investigation here will be more severely affected when a missing-p T cut is 
imposed. To take a specific example, the signal rate at benchmark point 1 becomes about 
48% of its original value when the requirement of a minimum missing-pr of 100 GeV is 
imposed. In contrast, at a a nearby point in the parameter space (mo = 200 GeV, 
m \/2 — 600 GeV, A = 100 GeV, tan/3 = 30, sgn(n) = +ve), the same missing-p^ 
cut allows about 97% survival of the jets + dimuons + jpx signal. Thus the response to 
missing-pr cuts turns out to be an effective tool of differentiation between our signal and 
that coming from MSSM, at least when R-parity is conserved. 
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3.2 Dimuon and two staus with two or more jets 

With the stringent demand on the hardness of the jets, this is a very clean signal, albeit less 
copious than the previous one. Such final states will require cascade decays of gluinos and 
squarks involving the charginos and heavier neutralinos. The same 'basic cuts' are imposed 
here, too, which are found to be sufficient in drastically reducing the SM backgrounds in the 
form of four muons together with two or more jets with > 100 GeV. 4 

In Figure 4 we plot the corresponding distributions for the above final states. It is the 
requirement of such energetic jets that removes the backgrounds to a large extent, as can 
be seen from Table 3. Furthermore, the plot of the scalar sum of the individual transverse 
momenta of the charged tracks in the muon chamber, in Figure 4(c), indicates that a cut of 
600 GeV on this sum washes out the backgrounds completely. The advantage of applying 
this cut is that there is no loss of signal events in the case of a relatively heavy NLSP. In 
Figure 4(a) and (b), the pr-distributions of the (harder) muon and the corresponding stau- 
track are seen to have a substantial overlap. Therefore, a distinction between them based 
on the thickness of the tracks as well as the information provided by measurement of the 
'time-of- flight' can be useful here. 



Final States 


Background 


Benchmark point 1 


Benchmark point 2 


2fi + 2fi 


83 


689 


103 


2fi + 2/x + (> 2) hard jets 


29 


686 


103 


2fi + 2fi + (> 2) hard jets 
(E Pt > 600 GeV) 





553 


89 



Table 3: The expected number of events for the signal and background with the different cuts 
imposed on the selection of events. Ept corresponds to the scalar sum of the individual transverse 
momenta of the charged tracks in the muon chamber. Choice of integrated luminosity is the same 
as Table 2. 



4 It should be noted that there are some uncertainties in the way PYTHIA models jet production, beyond 
two-jet final states. However, we retained PYTHIA results here, partly because the main point we wish to 
illustrate is unaffected, and partly due to the want of a completely reliable jet-production algorithm including 
SUSY. 
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4 Summary and conclusions 



We have studied a SUSY scenario with a stau-NLSP and an overwhelmingly right-chiral 
sneutrino as the LSP, where the sneutrino is at least partially responsible for the cold dark 
matter of the universe. A mass spectrum corresponding to such a scenario can be motivated 
in a SUGRA framework, including right-chiral sneutrinos whose masses remain practically 
frozen at the universal scalar mass at the SUSY breaking scale. We find that the superparticle 
cascades culminating into the production of stau-pairs give rise to very distinct signals of 
such a scenario. Although the charged tracks of the quasi-stable staus tend to fake muonic 
signals in the muon chamber, our analysis reveals considerable difference in their kinematic 
characters. Such difference can be used in a straightforward way to distinguish between the 
long-lived staus and the muons, and also to eliminate all standard model backgrounds. Since 
the mass spectrum under consideration here is as probable as one with a neutralino LSP in 
mSUGRA, further study of all possible ways of uncovering its signature at the LHC should 
be of paramount importance. 
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Figure 4: Kinematic distributions for the signal 2t\ + dimuon + (> 2) hard jets with basic cuts. In 
(a) and (b) the transverse momentum distributions for the harder f\ and harder muon of the signal 
is shown and (cKinematic distributions for the signal 2f\ + dimuon + (> 2) hard jets with basic 
cuts. In (a) and (b) the transverse momentum distributions for the harder f\ and harder muon 
of the signal is shown and (c) shows the scalar sum of px 's of ( dimuon + 2f\ ) of the signal. We 
follow the same notation as in Figure 2.) shows the scalar sum of px 's of ( dimuon + 2f\ ) of the 
signal. We follow the same notation as in Figure 2. 
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